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Abstract 



!> ' The non-singular, oscillating Friedman cosmology within the framework of 

t> ; 

General Relativity is considered. The general oscillatory solution given in terms 

m 

of elliptic functions and the conditions for its existence are discussed. It is shown 

in ' 

that the wall-like-matter and the small, but negative cosmological constant are 
required for oscillations. The oscillations can , in principle, be deep enough to 
allow standard hot universe processes like recombination and nucleosynthesis. 
^ ■ It is shown that the wall-like-matter and string-like-matter can be inter- 

c5 ' preted as scalar fields with some potentials. This may give another candidate 

for the dark matter which may be compatible with observational data. For an 
exact elementary oscillatory solution it is shown that the associated scalar field 
potential is oscillating as well. 
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1 INTRODUCTION 



The idea of a quasi-static and ever-lasting nonsingular universe seems still to be very 
attractive among cosmologists. It was first Einstein M who did not believe his equa- 
tions could really give in general non-static solutions, until the clear arguments of 
Friedman ||. With the discovery of Hubble redshifting and the mathematical work 
of Hawking and Penrose 0] about the inevitability of singularities rather few atten- 
tion was paid to non-singular models of the universe until the eighties. However, 
singularity-free Friedman models usually called bouncing models were considered 
early in thirties by Robertson ||. These models appear whenever the cosmologi- 
cal constant, introduced by Einstein to balance gravitational attraction, has a value 
between zero and a certain fixed number that gives the Einstein Static Universe 
< A < A st ||. The universe begins with some minimum value of the scale factor 
and then expands forever. Since the idea of inflation was established [[?], [J, the situ- 
ation changed because the scalar field responsible for the scenario necessarily had to 
break the energy conditions giving the effective negative pressure (repulsion) so much 
demanded by Einstein. 

If we raise the question what are the ways to avoid singularities one can distinguish 
a couple of approaches. In some opinions one should consider a quantum gravity as an 
appropriate gravity theory at high densities. However, most of the approaches refer 
to the alternative to General Relativity classical gravity theories. Worth mentioning 



is the early proposal made by Cartan |9], [10], |TTJ of nonsymmetric connection gravity 
theory with torsion in which both isotropic and anisotropic universes start with a very 
small but non-zero size |12|, [13| . Petry [|14], pl| considered a covariant gravity theory 



in which isotropic universes are always bouncing |IB . In the nonsymmetric metric 



1 



gravity theory of Moffat ]T7|, [R| |19j the black hole is replaced by a superdense object 



|20|| . Then, due to a repulsive force generated naturally in the theory, non-singular 
cosmological solutions should also be a rule. Other examples are the higher gravity 
theories such as: the second-order R + eR 2 theory ]2T| with a bounce solution for 
negative curvature |22|,|2§, the fourth-order R + XR^R^ / R + rR a()l& R a ^ s / R theory 



24j and Brandenberger's |25[ important proposal in which all isotropic solutions are 
the de Sitter, though non-singular, type. 

If we stay within General Relativity we can keep the simple Robertson- Walker 
geometry and either assume a specific equation of state for Friedman models |26|, 
or explore the scalar field coupled to gravity. Both approaches are actually equivalent 
since the scalar field can mimic different equations of state |28, p9| . 

The aim of this paper is to consider singularity-free bouncing Friedman universes 
within General Relativity admitting some fraction of negative pressure matter usually 
believed to describe cosmic strins and domain walls which we will be calling from now 
on, for the sake of some generality, string-like- matter and wall-like- matter The 
physical motivation for taking exotic matter into consideration is that we still do not 
have a solution of the dark non-baryonic matter problem. One can strongly believe, in 
the context of particle physics, that the exotic matter may be a very good candidate 
for the dark matter. If it appeared that the exotic matter existed one would believe 
in the non-singular universes of this paper. However, there exixts some objections 
referring to the compatibility of the domain walls with observations [31, |32| . Although 



some other suggestions emerged |33|, |34j, it is believed that even if walls formed in 
the early universe, they should decay by the present |35], |36fl . In order to avoid the 
domain walls and perhaps to give more freedom in compatibility with observations 
we replace the string-like-matter and wall-like-matter by scalar fields. We postulate 
that scalar fields may be the candidates for the dark matter in the present era of the 



universe f37| , Pq| . Our procedure is similar to the procedure applied for inflationary 
universes [[38] . First we assume the exotic equation of state and then we derive the 
form of the scalar field potentials. 

Among non-singular solutions we select those which are oscillating i.e. those 
for which the scale factor oscillates between certain fixed values. Such solutions 
first mentioned by Harrison [39] and recently discussed by Kardashev [32| originate 



from a qualitative change for bouncing models || after taking the negative pressure 
matter into account and admitting negative cosmological constant. These solutions 
have both contracting and expanding phase, so they may be compatible with the 
astronomical data. They last infinitely in time which might be satisfactory at least 
from philosophical point of view. 

Oscillating models might also have another advantage. Since they are quasi-static 
and require a similar balance to static models, then, according to the analysis of 



Gibbons [HJ [IT], they should posses very large entropy with value close a maximum 
admissible for the Einstein Static Universe. 

The plan of the paper is as follows. In Section 2 we discuss the existence of oscil- 
lating Friedman solutions and give the general exact oscillating solution in terms of 
the elliptic functions which has not been studied before. Also, we discuss a possibility 
to get very deep oscillations (i.e. those which can reach very high density to allow 
at least some standard hot universe processes) and to reduce the amount of exotic 
matter in form of wall-like-matter. In Section 3 we present the exact elementary non- 
oscillatory solutions which complete the discussion of Section 2. In Section 3 we give 
an alternative scalar field interpretation for exotic fluids necessary to force the uni- 
verse to oscillate. We choose a very simple model of the scalar fields of which potential 
and kinetic energy are proportional, following the discussion of Barrow and Saich [[42] . 
For a monotonic solution and the oscillating Kardashev's solution we present exact 



scalar fields and their potentials. It appears that for oscillating solutions for the scale 
factor the potentials are also periodic in the scalar fields. In Section 5 we comment 
on the results. 



2 OSCILLATING UNIVERSES 



Following our earlier discussion of exact analytic solutions of the Friedman equation 
||43| . fPfl we generalize it to the case when other negative pressure fluids are present, in 
particular wall-like-matter whose energy density scales as R~ l (R is the scale factor) 
|30| . The Friedman equation has the form 



^) = C r + C m R-k'R 2 + C w R 3 + ^R 4 , (2.1) 



where r is the conformal time defined by the cosmic time as 

dr = * (2.2) 

and k' = k — C s , k is the curvature index, and the constants C r , C m , C s , C w are 
constants responsible for the density of radiation, nonrelativistic matter, string-like- 
matter (whose energy density scales as R~ 2 ) and wall-like-matter respectively, and A 
is the cosmo logical constant. The Friedman equation (2.1) in terms of the so-called 
reduced blackbody temperature Q 

T = A^iT\ (2.3) 

where 

Ac h = ^C m , (2.4) 

2 The reduced blackbody temperature T(t) is proportional to the temperature of the microwave 
background if a ^ 0, and it is only interpreted as an inverse function of R(t) if a = [EHI 



is p| p 

(^j = a T* + ^-k'T 2 + [3T + ± (2.5) 
where the dimensionless parameter 

(3 = C w A c h (2.6) 

is responsible for the density of wall-like- matter (a = C r A c , A = A/A c ). In order to 
make qualitative analysis of the solutions of (2.5) we use the method of the associated 
mechanical system H^, In fact (2.5) can be considered as the energy equation 
of a one-dimensional mechanical system with T as a coordinate, the kinetic energy 
(dT/dr) 2 > 0, the potential 

Vv, a A T ) = - \ T * + k ' T2 -P T = -Qk>,a,f3,x(T) + (2-7) 



where 



Qk>,aAx(T) = «T 4 + ^T 3 - k'T 2 + /3T+^, (2.8) 



and the total energy |. 

The general shape of the potential (2.7) depends on the values of the factors 
k', a and f3. Instead of the full discussion of the solutions (2.5) we will concentrate 
on special cases which allow the universe to oscillate in time between the two fixed 



values of the scale factor R(t) (or T(r)) [32]. In order to have oscillations at least two 



kinds of matter are necessary i.e. the cosmological constant and the wall-like-matter. 
To make our considerations easier we put a = (no radiation) [] 



One can consider the full formula (2.7) that includes radiation as well, and the roots of the 

3 J 



cubic equation -aT 3 — |T 2 + k'T — (3 which we will call T\, T% and T3 have to fulfil the conditions 



Ti+T 2 + T 3 = ~, 
6a 

TiT 2 + TiTg + T 2 T 3 = — , 

a 
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The potential (2.7) is then 



--T 3 + k'T 2 - f3T } 
3 



(2.9) 



and it has one double root Tip — jk 1 for (3 = |& and two real roots for (3 < |fc as 
well as the root at T = 0. Its corresponding extrema are a double one at Tip = \k' 
for p = \k' 2 (which is an inflection point as well) and two real ones for [3 < \k' 2 . 
The general solution of (2.5) given in terms of the Weierstrass elliptic V function 



is 



T(t) 

where V is defined by the equation 
and the invariants 



12 



QaV(r) 



| — k'a 



(2.10) 



(TP 
dr 



AV 3 - g 2 V - g 3 , 



k' 2 aX 

q-2 = 1 , 

y 12 3 6 

,3 = 6 -3(^_ 2A _ 3 ^)-|(^ + f 



(2.11) 
(2.12) 



For a = case (no radiation) corresponding to the potential (2.9) we have the general 
solution of (2.5) in the elliptic form as well 



T{t) = 6 



(2.13) 



with the invariants given by (2.11)-(2.12) for a = and the discriminant [43 



Afc'AA = 92 - 27<? 3 2 = 2" 4 3 



4 Q -3 



4 

TiT 2 r 3 



P 2 (jk' 2 -2p)+k'\ (k l2 -?,p)-X 



(2.14) 



1 
a 



It follows from these relations that there should be at least one negative real root and the two real 
or complex conjugate roots, which means that radiation pressure can only change the left branch 
of the curve of the potential V at y t p(T) (i.e. negative T(t) - cf. Figs. 1-3) which is not physically 
relevant. 
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i.e. 

A fcW = -2- 4 3- 3 (A - A+) (A - A_) , (2.15) 
where the critical values of A are 

A ± = Ijfe' (k 12 - 3/3) ± ^k' 2 (k 12 - 3(3f + 4/3 2 Qfc' 2 - . (2.16) 

The associated mechanical system method is as follows. Having the exact shape of 
the potential (2.7) (Figs. 1-3) we cut the curves along the constant A lines. The point 
associated with the universe can only move in the upper part of the plane and we 
look for the possible solutions. We consider only the particular cases given by (2.13) 
where oscillations are possible (Ti j2 zeros and T lj2 extrema of (2.9)). These are: 

A. If k' > and < § we have 

Tifl = \{k' T \^k' 2 -2A(3), 
T 1)2 = \ (V =p sjk' 2 - 2/3) , 
and the oscillations corresponding to a constant A lines for 

A_ < A < 

with A_ given by (2.16) are possible in the central well of the potential (2.9) (Fig.l). 

B. If k' > and ^ = § we have 

7\ = T 2 = f 2 = ^k', 

Tl = T 

and the oscillations are possible for 

- Ik' 3 = A_ < A < A + = 0, 
8 
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i.e. in the central well of the appropriate potential (2.9) (Fig. 2). 
C. If k' > and § < ^ < \ we have 

H = T 2 = 0, 
1 



2i, 2 = -[k'T^/k' 2 -2(3 
and the oscillations are possible for 

A_ < A < A + < 0, 

with A-t given by (2.15) (Fig.3). 

There is apparently a double extremum for T — — if A — |, but it is also an 
inflection point of (2.9) and this case does not allow any oscillating solutions at all. 

In all considered cases A, B,C the discriminant (2.15) is positive and the funda- 
mental periodicity cell is a rectangle (Fig. 4) [j4"7fl . The roots e^, e 2 , e 3 of the equation 



43/" + tfjj/ + 03 = (2.17) 

are all real and since e\ + e 2 + e 3 = 0, at least one of them must be negative. One of 
the elementary periods uj is pure real and the other J is pure imaginary. The roots 
are given by the relations 

P(u) = e u 
V(oj + J) = e 2 , 
V{J) = e 3 , 

and 

e 1 <e 2 < e 3 . (2.18) 

However, since we consider the solutions of the equation for T(r) (cf. (2.5)) then we 
should define the roots of the equation 

2 A 

-T 3 — k'T 2 + (3T H — = 0, (2.19) 

3 3 



which are 

- 

T mm - 6e 3 + -, 
k' 

T m ax 6C2 ~\~ ~, 

k' 

Trecol 6C1 -|- —j 

and 

T m in ^ T m ax <C Trecol i (2.20) 

where T reco ; refers to a minimum of T(r), i.e. a maximum of R(t) for the recollapsing 
model associated with (2.13) and T m i n ,T max refer to a minimum and a maximum of 
T(t) for the oscillating model associated with (2.13). From (2.19) we conclude that 
T m in, Tmax and T reco i must be real and positive since T min + T ma x + T reco i = § and 
TminTmaxTrecoi = —A/2 and A is negative in the cases A, B, C. 

The general oscillatory solution for all the cases A, B, C can be expressed in terms 



of the Weierstrass £ function |45 



1 *I = -J- + \ T iC(r - r d ) - C(r + r d ) + 2((T d )} , (2.21) 



T(t) T m i n V A 

and the expression for the cosmic time, from (2.2), is 



3) i(T > 



<J / m.i.rj. 



+ ln^4, (2.22) 
tffa + r) 



so for r = t(r) = and T = T min (Fig.6). In both formulas (2.21) and (2.22) 
the conformal time is real but the zeros of the function T(r), namely — r d and are 
imaginary. Also, A in these expressions should be taken negative according to the 



general results of the existance of the oscillatory solutions and then y 3/A and J A/3 



9 



are imaginary. The periods are given by ||7| 



U! 



10 



dT 



Tmin J ~(T — T min )(T — T max ){T — T reco i) 

ecol dT 



T max (t _ T min ) (T - T max ) (T - T, 



recol I 



so uj is pure real and <J is pure imaginary. 

Let us now consider another type of the exact oscillating universes, the ones in 
which there is no dust (i.e. C m = in (2.1) and a = 0). Because of the definition 
(2.4) we cannot now use the relation (2.3) to work out the limit C m — > from (2.5). 
Instead we use the method by putting 

1 



M(r) 



R(r) 



so (2.1) becomes 



dr 



k'M 2 + C W M + 4- 



The potential analogous to (2.7) now is (Fig. 7) 

V kl , Cw = k'M 2 - C W M, 



(2.23) 



(2.24) 



(2.25) 



and provided k' > we have oscillations for A < 0. In this case the solution is 



elementary and oscillating [02 j. After a deparametrization based on (2.2) the exact 
solution can be written down as 

3 



R(t) = - 



2A 



AsmtJ-^ + C w 



where 



C v , > A. 



CI > 



A < 0, 



'C 2 w + -Ak' 



(2.26) 



(2.27) 



-Ak', 
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so the universe oscillates between R m in = — (3/2)A(— A + C w ) and 

Rmacc = -(3/2)A(A+C w ). 

From the shape of the associated potential (2.7) for (3 = |/c' 2 (Case B, Fig. 2) one 
can conclude that for ^ 

A_ = --k' 3 < A < 
8 

the universe will oscillate between the fixed extrema of R(t) and the smallest admis- 
sible value of the scale factor will be slightly larger than (cf. Section 3) 

Ru{r) = -^ = 2^, (2.28) 
AiT 2 k ' 



so the mass density of nonrelativistic matter will be [43 



3c 2 k' 3 3c 2 

Qm = ^C m R~ u 3 = cl^a- (2-29) 

If we assume k' = k — C s ~ 1 (only k = +1 is possible here, since k' > 0) which 
means that the amount of string-like-matter is small and most of the exotic matter 
is in the form of wall-like-matter , then from (2.28) 

g m ~ C m 2 ■ 2 ■ 10 26 ^. 

cm 

In order to calculate C m (with the dimension [cm]) we should apply the conservation 
law (2.28). For the maximum density at the nucleosynthesis scale 



g m ~ 10* 

cm 6 



we need 



C m ~ 1.4 ■ 10 n cm 



In fact the restrictions given below are actually the strongest ones since the minimum size of 
the oscillating universes in cases A, B, C is always bigger (i.e. T 2 < |fc' cf. Eq.(3.2)). 
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which can be achieved for instance either for the present nonrelativistic matter mass 
density g m o ~ 1.1-10~ 45 ^- (rather tiny) and the present radius Ro ~ 5-10 9 light years 
or for the present mass density g m o ~ 1.1 • lO -35 ^- (larger) and the smaller present 
radius Rq ~ 5 ■ 10 6 light years. However, for the maximum density at recombination 
scale 



we need 

C m ~ 1.4 • 10 22 cm, 

and it requires for example the present nonrelativistic matter density g m o ~ 1.1 • 
10~ 34 ^3 and the present radius R ~ 5 • 10 9 light years which is quite realistic. This 
means that oscillations might allow at least some standard early universe processes 
like nucleosynthesis or recombination. 

The amount of wall- like- matter can be reduced if we lessen the value of k' = k — C s 
(i.e. decrease the appropriate value of C m - cf.(2.28)), since the density of walls is 
proportional to k' 

P = \k' 2 = \{k-C s )\ 

and the wall-like-matter is replaced by string-like-matter ( < C s < 1 ). However, 
in such a case both the maximum and minimum of the potential (2.7) are approach- 
ing zero - the well becomes smaller which means that the appropriate periods of 
oscillations are smaller. 

Since the amplitude of oscillations corresponds to the length of the A = const, lines 
contained inside the well of the potential (2.7) (Figs. 1-3) one can easily notice a very 
interesting feature of the models A and B. It is that the closer to zero (necessarily 

12 



negative) value of the cosmological constant A, the bigger the maximum, and the 
smaller the minimum size of the universe can be, which means the longer period of 
oscillations. On the other hand, since we do not experience small period of oscillations 
(they can appear for large negative A - Figs. 1-3) we could obtain some restrictions on 
the possible value of the cosmological constant in these models. 

Finally, it is very interesting to mention that if we admit some fractional equations 
of state 

p= (e-l)p 

with s = ~n (n is non-integer) to derive the new Friedman equation different from 
(2.1), then a larger variety of possibilities for oscillations appear |p9, 48|. It can be 



explained in terms of the associated mechanical system simply because more wells 
of the potential are possible in which the universe can oscillate. Such fractional 
equations of state are possible since the exotic matter may have an effective equation 
of state anywhere in a range between well established values. For instance for strings 
| < e s < | and for walls \ < s w < I depending on their velocities 



3 NON-OSCILLATING EXACT UNIVERSES 



In this Section we concentrate on some particular solutions given by elementary func- 
tions. None of them is oscillatory but they complete the discussion of the solutions 
from the cases A, B, C as far as they can give some insight into the nature of the 
oscillatory solutions as well. 



5 It should be pointed out that string or wall stress is in fact anisotropic and might be considered 
by using an anisotropic model of spacetime rather than the isotropic one [ |50| ]. However, it seems that 
there is no reason for strings or walls to be formed anisotropically throughout the whole universe 



and we can suitably average the network over all directions |51 
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Generally, if the discriminant (2.14) vanishes the solution is elementary. From 

(2.15) one can realize that it happens for A = A + and A = A_ with A± given by 

(2.16) . With the condition of the vanishing discriminant the equation (2.5) factorizes 
to the form 



with 



T = T 1>2 = \{k'T \A' 2 - W) (3-2) 



with A = const. With given by (3.2) we have two Einstein Static Universes in 
the cases A and C (Figs.l and 3). One of them should be stable because it lies at the 
bottom of the well of the potential (2.7). In the case B i.e. for — | (3.2) reduces 
to 

T = T 1>2 = % ( 1 =F I) ■ (3.3) 



2 V 2, 

With every stable solution for T\ it is also associated the elementary recollapsing 
solution for A = A_. On the other hand with each unstable solution for T 2 there are 
also associated two asymptotic solutions for A = A_. 

Since the types of the elementary solutions of (3.1) in the cases A and C are 
analogous to the solutions in the case B with only numerical difference between (3.2) 
and (3.3) we will be discussing the latter case which is mathematically simpler. The 
stable Einstein Static Universe appears for A = A_ = — |A;' 3 and 

Rs = 4^ = 6%. (3.4) 
A|Ti k 

The equation (3.1) for A = A_ reads as 

(£HH)V-*% 
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and its solution for T > k' describes the evolution of the closed universe from the 
big-bang to the big-crunch in the form 



R(r) 



MT(r) k! 



Ac' 2 

1 + ftan^vfr 



(3.6) 



The unstable Einstein Static Universe appears for A = A + = and 



R 



. Cn 



U 1 _ 

AIT 2 



The equation (3.1) for A = A + = reads as 



(3.7) 



'dT\ 2 1, 



— = -T T k 

dr J 3 V 4 



(3.8) 



There are two asymptotic to (3.7) solutions of (3.8). For < T < -k! we have 



R(r) 



MT(t) 



4 A -i , 2 1 k' 
—Ac 2 coth 2 -W— r, 
3k' 2V 2 ' 



(3.9) 



and according to (2.2) 



t(r) 



8V2 



(k ,3 A c y 



1 k' 1 k' 
-\ — r — coth -\ — r 

2 V 2 2 V 2 



(3.10) 



Then for T > \k' we have 



R(r) 



MT(r) 



4 — ^ l2 l fc' 



(3.11) 



and 



tr) 



8V2 



(k' 3 A, 



1 k' 1 k' 
- \ / — r — tanh - -t / — r 

2 V 2 2 V 2 



(3.12) 
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4 SCALAR FIELD COSMOLOGIES 



In this section we try to give an alternative interpretation for exotic matter appearing 
in the Friedman equation (2.1) in terms of the scalar fields and appropriate potentials. 
The motivation is that we still are looking for the new candidates for the dark matter 
in cosmology and we may think about scalar fields existing at the present era of the 
evolution as good candidates. Such an interpretation seems to give more freedom in 
imposing the observational constraints than string-like-matter and wall-like-matter 



1 37J. Our procedure is a reverse of the standard procedure dealing with scalar fields. 
Instead of starting with an explicit potential to derive the equation of state, we start 
with the equations of state for string-like-matter and wall-like-matter (cf. Section 2, 
formula (2.1)) in order to derive the exact potentials ||38|| . The energy density and 
pressure for two mutually noninteracting scalar fields minimally coupled to gravity 
can be written as 

Pi=\$ + Vi(<pi), (4.1) 

Pi = ~$-VM), (4.2) 

where i = 1,2. If, apart from this, we admit a barotropic fluid with the equation of 
state p = [y — l)p we will obtain the rewritten Friedman equation (2.1) in the form 



3 (f ) =3H* = p+ I ft + I ft + V^) + V 2 {<p 2 ) - A (4.3) 
and the wave equations for both fields 

<p x + 3#0x + VIM = 0, (4.4) 
Cp 2 + 3H<p 2 + V 2 '(ip 2 ) = 0, (4.5) 
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where (. . .)' = and (. . .) 



■g-. Following Barrow and Saich (1993) we assume a 



simple hypothesis that kinetic and potential energies of the (pi,(f2 fields are propor- 
tional, so 

1 



ttiVi(y?i) = -cp v 
"2^2(^2) = ^l, 



(4.6) 
(4.7) 



with a±, «2 = const.. If we use the conformal time (2.2) then we can write down the 
solutions of (4.4)-(4.5) under the conditions (4.6)-(4.7) in the following way 



(p i>T = AiR 



(4- 



where Aj are constants of integration. The rewritten Friedman equation (2.1) becomes 
(j = 1,2,3) 



'dRY 



CjR-^+ 4 + - (l + —) A 2 R^ - kR 2 . (4.9) 

2 V oij J 



Comparing (2.1) with (4.9) we realize that in order to exchange string-like-matter 
and wall-like-matter for scalar fields we have to take |H4 



°> C ' = J> 

1, O2 = C m 
4 

- , C 3 = C r 



ocx = \,C s =*-A\^A x = J 2 ~{k>-k), 



(4.10) 



a 2 = -,C tt = 3A!=*A 5 
5 



\ 



T 



(4.11) 



In order to obtain the asymptotic solution (3.9) the exact scalar fields should be 



tpi = Air + cp 01 
17 



(4.12) 



for a\ — | and 



Lf 2 = In 



V?02 



cosh —\ — r 
2V 2 



2A- 



(4.13) 



for c*2 = |, with ^01,^02 = const, and Ai, A 2 given by (4.10)-(4.11). We can thus 
write down the required solutions for the scalar fields as 



<Pi 



V?2 



- (k' - k)r + <p 01 , 



In 



V?02 



cosher 



(4.14) 



(4.15) 



with the associated potentials 

A 2 

K(r) - % 

V,(r) = ^ 



^(fc'-A;)A;' 2 A c coth 2 iJ|r, 



15 | , 1 /fc' 

=A£ coth -W —r, 

2 i? 16^2 2 V 2 



(4.16) 



(4.17) 



with cosmic time t(r) given by (3.10). The potentials as functions of the scalar fields 
are (Fig. 8) 



V 2 {<P2) 



- (k 1 - k) k' 2 A c coth 



J ) A;' 2 a! coth 



16 



cosh 1 (p 2 exp 



' 9ip 2 
A\fk' 



(4.18) 



(4.19) 



Of course one could also obtain exact potentials for exact models (3.6) and (3.11), 
but they do not oscillate as well. The general solution for oscillating models is given 
in terms of elliptic functions by (2.10) ( or (2.13) if a = 0) and the exact potentials 
could be calculated using Eq.(4.8). As another important example we consider the 
elementary oscillating model given by p2 |. 



According to (2.2) the Friedman equation (4.9) in terms of cosmic time t(r) is 



'dR\ 



dt J 



\ = C 3 R~^ +2 + - f 1 + - ) AjR 



1 



( 2 ^) - 



k. 



(4.20) 



The exact solution given by |E| for the case without matter and radiation is given by 
(2.26), the values of constants a>\ and a>2 are given by (4.10)-(4.11) and the resulting 
fields are 



2 / k\ 
a - 1 — — arctan 

\ 3 V k J 



■lAk' 



(4.21) 



but it reduces to a constant, if there is no strings (i.e. if k' — k — +1), and 



2 a a, 



dt 



C w + A sin t\ — 



tan 



A 

T 



, (4.22) 



where V 1 is the inverse function to the Weierstrass elliptic V function with invariants 
given by 

92 = \(ci-AAk' 



1 ( C 2 
93 ~ -g 6 - I "g- 



+ 4AA;' J 



(4.23) 
(4.24) 



It is easy to check that (4.21) and (4.22) are oscillatory in t. According to (4.6), (4.7), 
(4.10) and (4.11) 



Vi(t) 



v 2 (t) 



2 k 



\3\ k'jAVtt 



2A 2 fa, 



cos 2 \t 



C w tan 



+ A 



2 „ ' 
Uk' 



(4.25) 



3 V 3 c w + AsmtJ-± 



(4.26) 



From (4.25) and (4.26) we realize that the potentials are oscillatory in t i.e. V\ 
oscillates between B(-A/?,kk')(C w - A) and B(-4/3Ak')(C w + A), where 
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B = 4- 1 A;'-2,/2/3(l - k/k') and V 2 oscillates between 2A 2 /3JC w /3/(C w - A) and 



2A /3\/C w /3/(C w + A). If we express explicitly potentials in terms of fields we have 



V 2 (ip 2 



k'~ 2 Cr 2 



32 J\ (k' - k) 



C 2 W + 



Ak' tan 2 $'_ k) {(fi - (poi) 



-A 



l C w 1 


3 C w + A sin {2 arctan 




]} 



(4.27) 



(4.28) 



The potential Vi(ipi) oscillates between Vd = — |Afc' and V u = 4C^+4Afc'. A prove 
that ^2(^2) is also oscillating in (p 2 is more complicated. The discussion involves the 
analysis of the discriminant A = — g% with g 2 and #3 given by (4.23)-(4.24) and 
then the analysis of the solutions similar to that given in Section 2 for the general 
oscillatory solution subject to the conditions (2.27). Finally, the result is that Vi(ipi) 
and V 2 (ip 2 ) in (4.27)-(4.28) do oscillate in cpi and <p 2 . 

These results together with (4.18) and (4.19) may suggest that at least under the 
assumption of proportionality of kinetic and potential energies (4. 6)- (4. 7) the mono- 
tonic type of the solution for the scale factor R(t) leads to the monotonic dependence 
of the potential V on the scalar field (p and the oscillating type of solution for R(t) 
leads to the oscillating dependence of the potential V on the scalar field (p. 



5 DISCUSSION 



In this paper we discussed a class of nonsingular oscillatory universes following Kar- 



dashev's [S2| discussion of admitting essential fraction of stable domain walls as the 
exotic matter and the negative cosmological constant. In the context of the dark non- 
baryonic matter problem both kinds of fluids can serve as dark matter candidates. 

20 



The former acts as a source of repulsive gravity but the latter acts as a source of 
attraction. The balance between these two sources allows the universe to fall into 
oscillations. In fact, this is very similar situation to the case of the Einstein Static 
Universe, where there is a balance between nonrelativistic positive pressure matter 
and repulsive positive cosmological constant. According to our analysis the Einstein 
static models exist as a result of a balance between the exotic matter and the at- 
tractive negative cosmological constant. Unlike the original Einstein Static Universe 
these models seem to be mechanically stable as they lie at the bottom of the asociated 
mechanical potential (2.7), but their stability should be examined more carefully. 

Since the value of the negative cosmological constant can be reduced to be very 
close to zero for oscillatory solutions (cf. Section 2), the most severe problem of 
possible compatibility with observations refers to exotic matter, especially to domain 



walls PHI . Because of that we proposed a replacement of the exotic matter by scalar 
fields which may serve as dark matter candidates as well. We inversed the standard 
procedure - first we had the equation of state and then we derived exact potentials 
P%| . The potentials we studied seem to keep the same properties as the solutions for 
the scale factors. It means that monotonic solutions for R(t) (cf.(3.9)) give mono- 
tonic potentials (cf.(4.18)-(4.19)) and oscillatory solutions (cf.(2.26)) give oscillatory 
potentials (cf.(4.27)-(4.28)). 

It is worth emphasizing that the paper is based on many simplifications and the 
question is what are the properties of oscillatory solutions in more physically realistic 
situations. Firstly, in the discussion of oscillatory solutions of Section 2 we have 
neither considered any non-isotropic geometry nor non-isotropic fluid description for 
instance for string-like-matter and wall-like-matter. Secondly, in the scalar fields 
interpretation of Section 4 we have assumed a very simple proportionality relation 
for the kinetic and potential energies, the minimal coupling of the fields to gravity 
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and the lack of interaction between them. Of course all these assumptions might 
not necessarily apply and the validity of the presented results should be considered 
towards these points. 

Especially, in future work one can consider some non-isotropic cosmologies Q 
and non-minimally coupled fields |55], but these approaches will certainly be 
much more difficult to treat analitycally. Also, one can discuss oscillating quasi-static 



universes in the context of the maximum entropy analysis given by Gibbons [}40| , [41 
Finally, we should mention one common feature of our results with the results of Petry 
|i~6|1 in his conformal gravity theory. Namely, in both theories nonsingular oscillating 



solutions appear for the negative cosmological constant. How far this results could 
be true for other gravity theories is a matter for further consideration. 

Acknowledgments 

The author wishes to thank to John Barrow, Arne Larsen, Jerzy Stelmach and 
David Wands for careful reading of the manuscript and helpful suggestions. 



22 



References 

[1] A. Einstein, Sitzber. Preuss. Akad. Wiss. 1 (1917) 142 

[2] J. D. Barrow k P. Saich, Class. Quantum Grav. 10 (1993) 279 

[3] A. Friedman, Z. Physik 10 (1922) 377 

[4] S. W. Hawking & G. F. R. Ellis (1973) The Large-Scale Structure of Spacetime 
Cambridge University Press 

[5] H. P. Robertson, Rev. Mod. Phys. 5 (1933) 51 

[6] M. Rowan-Robinson, (1981) Cosmology Oxford University Press p. 121 

[7] A. H. Guth, Phys. Rev. D23 (1981) 347 

[8] A. D. Linde, Phys. Lett. B108 (1982) 389 

[9] E. Cartan, C. R. Acad. Sci. (Paris) 174 (1922) 593 

[10] E. Cartan, Ann. Ec. Norm. Sup. 40 (1923) 325 

[11] E. Cartan, Ann. Ec. Norm. Sup. 41 (1924) 1 

[12] A. Trautman, Nature Phys. Sci. 242 (1973) 7 

[13] W. Kopczyhski, Phys. Lett. A43 (1973) 63 

[14] W. Petry, Gen. Rel. Grav. 10 (1979) 599 

[15] W. Petry, Gen. Rel. Grav. 13 (1981) 865 

[16] W. Petry, Gen. Rel. Grav. 13 1981) 1057 

[17] J. W. Moffat, Phys. Rev. D19 (1979) 3554 

23 



[18] J. W. Moffat, Phys. Rev. D19 (1979) 3562 

[19] J. W. Moffat, Journ. Math. Phys. 21 (1980) 1978 

[20] J. W. Moffat, University of Toronto Preprint UTPT-84-14 

[21] J. D. Barrow & A. C. Otewill, J. Phys. A16 (1983) 2757 

[22] D. N. Page, Phys. Rev. D36 (1987) 1607 

[23] D. H. Coule & M. S. Madsen, Phys. Lett. B226 (1989) 31 

[24] C. Xu & G. F. R. Ellis, Found. Phys. Lett. 5 (1992) 365 

[25] R. Brandenberger, V. Mukhanov & A. Sornberger, Phys. Rev. D48 (1993) 1629 

[26] M. Israelit k N. Rosen, Astroph. Journ. 342 (1989) 627 

[27] M. Israelit & N. Rosen, Astr. Space Sci. 204 (1993) 317 

[28] S. §. Bayin, F. I. Cooperstook & V. Faraoni, Astroph. Journ. 428 (1994) 439 

[29] A. Liddle, Phys. Lett. B220 (1989) 502 

[30] J. Stelmach, R. Byrka & M. P. Dabrowski, Nucl. Phys. B406 (1993) 471 

[31] N. Turok, Phys. Rev. Lett. 63 (1989) 2625 

[32] N. S. Kardashev, Mon. Not. R. astr. Soc. 243 (1990) 252 

[33] B. G. Gelmini, M. Gleiser & E. W. Kolb, Phys. Rev. D39 (1989) 1558 

[34] G. Goetz, Journ. Math. Phys. 31 (1990) 2683 

[35] B. Rai & G. Senjanovic, Phys. Rev. D49 (1994) 2729 

[36] M. Cvetic, S. Griffies & H. A. Soleng, Phys. Rev. D48 (1993) 2613 

24 



[37] J. McDonald, Phys. Rev. D48 (1993) 2462 

[38] J. D. Barrow, Phys. Lett. B325 (1990) 40 

[39] E. R. Harrison, Mon. Not. R. astr. Soc. 137 (1967) 69 

[40] G. W. Gibbons, Nucl. Phys. B292 (1987) 784 

[41] G. W. Gibbons, Nucl. Phys. B310 (1988) 636 

[42] J. D. Barrow & P. Saich, Class. Quantum Grav. 10 (1993) 279 

[43] M. P. Dabrowski & J. Stelmach, Ann. Phys. 166 (1986) 422 

[44] M. P. Dabrowski & J. Stelmach, Astron. Journ. 97 (1989) 978 

[45] R. Coquereaux R & A. Grossmann, Ann. Phys. 143 (1982) 296 

[46] S. Chandrasekhar (1983) Mathematical Theory of Black Holes Oxford University 
Press 

[47] F. Tricomi (1937) Funzioni Ellittische Nicola Zanichelli Editore, Bologna 

[48] M. P. Dabrowski (1986) unpublished 

[49] E. W. Kolb & M. S. Turner (1990) The Early Universe Addison- Wesley 

[50] J. D. Barrow, G. J. Galloway & F. J. Tipler, Mon. Not. R. astr. Soc. 223 (1986) 
835 

[51] M. Hindmarsh (1994) private communications 

[52] N. Weiss, Phys. Rev. D39 (1989) 1517 

[53] A. Feinstein & I. Ibanez, Class. Quantum Grav. 10 (1993) 93 

25 



[54] A. Feinstein & I. Ibanez, Class. Quantum Grav. 10 (1993) L227 

[55] T. Futamase & K. Maeda, Phys. Rev. D39 (1989) 399 

[56] T. Futamase, T. Rothman & R. Matzner, Phys. Rev. D39 (1989) 405 



26 



Figure Captions 



Fig.l. The potential of the associated mechanical system (2.9) for the Friedman equa- 
tion (2.1) with k! > and (3 < |A;' 2 . The universe can oscillate in the central well 
between T = and T = 7\ = f (k' - ^jk' 2 - p). 

Fig. 2. The potential of the associated mechanical system (2.9) for the Freidman 
equation (2.1) with k! > and (3 = |A;' 2 . The universe can oscillate in the central well 
between T = and T = T\ = T 2 = f ' 2 = \k' . 

Fig. 3. The potential of the associated mechanical system (2.9) for the Friedman 
equation (2.1) with k! > and |/c' 2 < (3 < \k' 2 . The universe can oscillate in the 
central well between T = and T = T 2 = \(k' + ^k r2 - 2/?). 

Fig. 4. The elementary periodicity cell for oscillatory solutions. The period uo is 
pure real while the period J is pure imaginary. 

Fig. 5. The exact oscillatory solution in the cases A, B, C from Section 2. This 
solution is given in terms of the Weierstrass ( functions by (2.21). 

Fig. 6. The exact recollapsing solution that appears in the cases A, B, C of Section 
2 (formula (2.13)). 

Fig. 7. The potential of associated mechanical system (2.22) for Kardashev's case 
(no nonrelativistic matter and radiation). Oscillations appear for A < and k! > 0. 
Fig. 8. The exact potentials for the asymptotic solution (3.9). 
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